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PAOLO LELLA 

Abstract. Borel-fixed ideals play a key role in the study of Hilbert schemes. Indeed each component and each 
intersection of components of a Hilbert scheme contains at least one Borel-fixed point, i.e. a point corresponding 
to a subscheme defined by a Borel-fixed ideal. Moreover Borel-fixed ideals have good combinatorial properties, 
which make them very interesting in an algorithmic perspective. In this paper, we propose an implementation 
of the algorithm computing all the saturated Borel-fixed ideals with number of variables and Hilbert polynomial 
assigned, introduced from a theoretical point of view in the paper "Segment ideals and Hilbert schemes of 
points". Discrete Mathematics 311 (2011). 



1. Introduction 

The Hilbert scheme Hilbp^^-, parametrizes all the subschemes and all the families of subschemes of the 
projective space P" with Hilbert polynomial p{t). Borel-fixed ideals are a basic tool for the direct study of 
Hilbert schemes, because 

• each component and each intersection of components of a Hilbert scheme contains at least one point 
defined by a Borel-fixed ideals; 

• the Borel-fixed ideals have a strong combinatorial property which makes them very convenient also from 
an algorithmic perspective. 

For instance, Hartshorne in his thesis [ ] proved the connectedness of Hilbert scheme constructing sequences 
of deformations of Borel-fixed ideals (he called them balanced ideals) which lead from any point of the Hilbert 
scheme to the point determined by the unique saturated lexicographic associated to the Hilbert polynomial p(t). 

More recently, many papers [2, 3, 4, 5, 6, 10, 11] by Bertone, Cioffi, Marinari, Roggero and the author 
showed how to use Borel-fixed ideals for a local study of the Hilbert scheme, mainly constructing families of 
ideals sharing with a fixed Borel-fixed ideals the same basis of the quotient space. 

Therefore it was very important to have an algorithm computing for each Hilbert scheme Hilbp^^^ all the 
points defined by Borel-fixed ideals. An algorithm for computing all the Borel-fixed ideals in K[xo, . . . , Xn] with 
Hilbert polynomial p(t) based on a combinatorial approach has been proposed in ["] and another algorithm with 
slight differences to this, and published later, is presented in [12]. 

In this paper, we describe a concrete implementation of the algorithm, that turns out to be very efficient as 
we will show with an experimental analysis. 



2. Notation and general setting 

We will consider a field K of characteristic and for any polynomial ring K[a;] := K[a;o, ■ • ■ , Xn] we will order 
the variables as 2;„ > ... > xq. Following the notation of [ ], we will refer to the Hilbert polynomial of a 
homogeneous ideal / C K[a;o, . . . , x„] as the Hilbert polynomial of the graded module K.[xo, . . . , Xn]/I, i-e. 

p{t) =dimTss,K[xo,...,Xn]t/It, i > 0. 

A homogeneous ideal / is said Borel-fixed if it is fixed by the action of the Borel subgroup of upper triangular 
matrices. Looking at the action of the elements of the Borel subgroup Id + Eij, where i < j and Eij is a 
matrix with all entries equal to except for the entry of the i-th row and j-h column equal to 1, it is possible 
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to prove that the ideal / has a nice combinatorial property: in fact an ideal / is Borel-fixed if and only if it is 
a monomial ideal and 

As done in [5, 10], we define the elementary Borcl moves: 

• e~ as the element ^^^^-^ in the field of fraction K(x) of K[x]; 

• as the element 2i±i g K(.t); 

and for any monomial x" we will say that e~ (resp. e~^) is admissible on x" if €~{x°') = 2^i=J-x" G IK[x] (resp. 
e^(a;") = ^^+^ a;" S K[a;]). We will use additive notation to denote the composition of an elementary move with 

itself, for instance 2e7 = e7 o e7 = fii^^ . 

J J J \ Xj J 

■p(n, r) will denote the poset of all monomials in ]K[x]r with the Borel partial order <.b given by the transitive 
closure of the order relations 

e+(x") >BX" >se-(x") 

V x" G IK[x]r, V e^,ej admissible, and we call Borel set any subset of a poset V{n,r) closed under increasing 
elementary moves. By definition, given a Borel-fixed ideal / C K.[x], the monomial basis of each homogeneous 
piece /( of fixed degree will define a Borel set, thus we will write {If} meaning the Borel set of P{n, t) defined by 
the monomials in If. On the other hand, given a Borel set ^ C V{n^r), we will denote by {iiS) the Borel-fixed 
ideal generated by the monomial of ^ and by {ii§Y^'^ its saturation. 
For any monomial we denote 

• maxx" — max{xi s.t. Xi \ x'^}; 

• minx" = min{xj s.t. Xj \ a;"}. 

For the order chosen on the variables, it will not be misleading to write only the index of the variable instead 
of the variable itself. 

For any Borel set !3§ C 'P{n, r) we will denote by the complement set Tin, r) \ SS. Obviously SS'~' is closed 
under decreasing elementary moves and we will call such a set order set^ being its dehomogeneization (imposing 
Xq = 1) an order ideal. 

For any subset S of 'P{n, r), S(^^i) will denote the subset of S of the monomials with minimum greater than 
or equal to i: 

S(^i) :— {x" G S I minx" ^ i} . 

Obviously S'f^o) — 

As well known, given any Borel-fixed ideal /, its saturation is generated by the monomials obtained 
substituting = 1 in the monomials that generate /. Moreover, xq is not a nihilpotent element for any 
Borel-fixed ideals. 

Proposition 2.1. Let I C K[x] be a Borel-fixed ideal. The linear form xq is regular for I. Thus, for any shift 
of degree t, there is the short exact sequence induced by the multiplication by Xg 

^ - 1) ^ ffl(i) ^ ^ 0. (1) 

This regular sequence says that the ideal (/, xq) has Hilbert polynomial Ap{t) — p{t)—p{t—l). We will perform 
hyperplane section repeatedly so we define A^p{t) := p{t) and recursively A''p{t) :— A'^^^p{t) — A''^^p{t — 1). 
It is very easy to manipulate Hilbert polynomial considering the Gotzmann decomposition. 

Definition 2.2. An admissible Hilbert polynomial p{t) has a unique Gotzmann decomposition 
ai^ . . Or- The number r of terms in this sum is said Gotzmann number of p{t). 



It can be easily proved that the Gotzmann decomposition of A.p{t) can be obtained from the decomposition 
of p(t) decreasing by 1 each index ai and discarding the binomial coefficients with a negative index below. 
Moreover we define an inverse operator E that associates to p{t) the Hilbert polynomial obtained increasing by 
1 all the indices a^. 

Example 2.1. The Hilbert polynomial p{t) = M + 1 has Gotzmann decomposition 



Hence 



We conclude the background materials with the properties linking the Hilbert polynomial to the regularity 
of an ideal. 

Definition 2.3. A coherent sheaf J- over P" is said m-regular ij for every i > 

^{Tim-i)) =0. (3) 

The Castelnuovo-Mumford regularity of T is the smallest integer m, for which J- is m-regular. 

Theorem 2.4 (Gotzmann's Regularity Theorem). Let A be any K-algebra and let Z C Proj A[xo, . . . , a;„] 
be any subscheme with Hilbert polynomial p{t), whose Gotzmann number is r. Then the sheaf of ideals Xz is 
r-regular. 

Proposition 2.5 ([ , Proposition 2.6]). Let I be a saturated ideal in 'K.[x\. The regularity of L is equal to the 
Castelnuovo-Mumford regularity of the sheaf of ideals X obtained from the sheafification of L . 

Proposition 2.6 ([ , ]). The regularity of a Borel-fixed ideal L C K[a;o, . . . is equal to the maximal degree 
of one of its generators. 

3. Main properties 

Now we recall the technical properties that we need to guarantee the correctness of the algorithm, which are 
already proved in [ ] and that we here adapt at the current notation. 

Lemma 3.1. Let J C K[x] be a saturated Borel-fixed with Hilbert polynomial p(t) whose Gotzmann number 
is r. Let x^ be a minimal monomial for <b of {Jr} C V{n,r) such that mina;^ = Xq. Then the ideal 
I = {{Jr} \ {x^}Y'^^ is Borel-fixed with Hilbert polynomial p{t) — p{t) + 1. 

Proof. First, note that by definition of minimal monomial, {Ir} is still a Borel set. Called q{t) the volume 
polynomial of J, we show that / has volume polynomial q{t) — q{t) — 1 applying Gotzmann's Persistence 
Theorem [. Theorem 3.8], i.e. proving that 

dimK Jr - dimK Ir = dimx Jr+i - dimx Ir+i = 1- 

By construction dimjK Jr — dimK Ir = 1- The Borel condition ensures that x^xq G Jr+i \ Ir+i and there are 
no other elements, because x^xq is the only monomial that cannot be generated from the monomials in Ir by 
multiplication of a single variable. In fact let us consider the monomial XiX^ ^ i > 0. Since xq \ x^ the following 
identity holds: 

Xix'^ — • . . . • — xqx^ = e,_i o ■ • ■ o en (x'^ )xa 

X;-l Xq 

and for each j, ° ■ ■ ■ ° ^ti^^) belongs to Ir, by the minimality of x^ . □ 
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Proposition 3.2 ([ , Proposition 5.2]). Consider a saturated Borel-fixed ideal J C K.[x] defining a subscheme 
with Hilbert polynomial p{t) whose Gotzmann number is r. Let I ~ J|xi=2;o=i be its double saturation and let 
p{t) be the Hilbert polynomial of the subscheme defined by I in P". Then 

p{t) = p{t) - diniK It + diniK Jr- (4) 

Theorem 3.3. Let p(t) be an admissible Hilbert polynomial in P" and let r be an integer greater than or equal 
to the Gotzmann number ofp{t). There is a bijective function between the setX^^^s^ of the saturated Borel-fixed 
ideals of ]K[a;o, . . . , Xn\ with Hilbert polynomial p(t) and the set 



nn 



C 'P{n,r) Borel set s.t. 
^A'p{r), Vi = 0, 



More precisely 



'Tn ,1-1. KJn 

Pit) ^ ^ ^p(t) 

J {Jr} 

l^^Y^t ^ gg_ 

Proof. First of all note that if the two maps are well-defined, i.e. for each J, J^r — (Jr) (Proposition 2.6 and 
Gotzmann's Regularity Theorem), 

J {Jr} ({J.})''' - J, 

— > {ssY''^ — > {{s§)f^]=s§. 

Let J C K[a;] be a Borel-fixed with Hilbert polynomial p{t) and let ,yV = V{n,r) \ {Jr}- Obviously 
\^{^o)\ = \^\ — pi''') — A"p(r). Using the short exact sequence (f), we determine the Borel ideal / = 
(J, xqY^''^ n IK[a;i, . . . ,x„] with Hilbert polynomial Ap{t). Thus being {Lr} = {Jr}{^i) C V{n - I, r), j-^jti) | = 
— Ap{r). Since / is Borel-fixed in the polynomial ring ]K[a;i, . . . ,x„] we can repeat the reasoning with 
the hyperplane section defined by xi = and so on. 

Let us now consider a Borel set £§ C V{n,r), such that the complement yy = ^ satisfies the condition 
l^^i) I = A*p(r) for every i. Firstly reg {{^Y'^^) ^ r by Proposition 2.6, so let us prove that {S§Y^^ has Hilbert 
polynomial p{t). We proceed by induction on the degree d of the Hilbert polynomial. For any n, if degp(i) = 0, 
then yyi^i) = 0, for every i ^ 1, since Ap(t) = 0, that is for any x^ G o/K, minx^ = 0. Applying repeatedly 
Lemma 3.1 starting from the Hilbert polynomial p(t) = (corresponding to the ideal (1)), we obtain that {SSY^"^ 
defines a module K[2:]/ {!3§Y^'^ having constant Hilbert polynomial p{t){= p{r)). 

Let us know suppose that the map S§ — iJ^Y^"^ is well-defined for any Hilbert polynomial of degree d — 1 
and let p{t) be a Hilbert polynomial of degree d. SS — ^(-^i) C ■P(n — l,r) realizes the condition of the 
theorem w.r.t. the Hilbert polynomial Ap{t) and degAp{t) = d — 1. Hence by the inductive hypothesis the 
ideal {^Y""^ C K[xi, . . . ,Xn] has Hilbert polynomial Ap{t). Let p{t) be the Hilbert polynomial of {^Y'^*' 
K[xo, . . . ,Xn]: p{t) — p{t) + a, because Ap{t) ~ Ap{t). {SSY^"^ turns out to be the xi-saturation of {S§Y^^ ^ 
so by Proposition 3.2 the Hilbert polynomial of ]K[a;o, . . . , a;„]/ {SSY'^^ differs by a constant from p{t) and since 
\.vV\ ~ |o/f(^o)l — p{t) it coincides with p{t). □ 



4. The algorithm 

Therefore to compute the saturated Borel-fixed ideals we can construct Borel sets with the prescribed prop- 
erty. The proof of Theorem 3.3 suggests to use a recursive algorithm: i.e. to determine the Borel sets in V{n, r) 
corresponding to the Hilbert polynomial p{t), we begin computing the Borel sets in Viji — l,r) corresponding 
to the Hilbert polynomial Ap{t). 

Let ^ C V{n — l,r) a Borel set corresponding to the Hilbert polynomial Ap{t) and let ,yV = 3§ . In order 
for ^ to be the restriction ^(^i) of a Borel set 3§ C 7^(n, r) (where 7^(n, r) contains one more variable smaller 
than variables in P{n — l,r)), each monomial that can be obtained by decreasing moves from a monomial in 
c/K has to belong to — SS'~ . This extension of an order set jV C V(n — 1, r) to an order set jV C P(n, r) 
has an ideal interpretation. 



Lemma 4.1 ([ , Lemma 5.1]). Let ^ C P(n — 1, r) be a Borel set and let JV — . Moreover let jV C V{n, r) 
he the order set containing the monomials in ,yV and all those obtained by decreasing moves from them. Then, 

.yV = V{n, r) \ [{{Wr' ■ K[xo, . • ■ ,x„]) J. (5) 

Proof. Let us call ^ the Borel set containing the monomials of degree r belonging to the ideal {^Y^*- ■ 
K.[xot ■ . ,Xn\. Let — a;"" ■ ■ ■ Xq° be a monomial of V{n,r) and suppose mina;" = 0, i.e. ao > 0. The 
monomial aoGo (a;") = a;^" •■• x"^^"" belongs to V{n—l,r), so either aoeJ(a;") S ^ or aoSQ (z") e ._yf . 

If x^" ■■■x°^+"° €W, then x^^ ■ ■ ■ x":^" is in {W)'"'^ and so x" £ {W^"^ -KIxq, . ■ . ,Xn], otherwise 
x^" • • • xi ^+"° e implies a;« G ^. □ 

By Proposition 3.2, we know that the Hilbert polynomial corresponding to a Borel set ^ of the type {{{^)^^^ ■ 
K[x])r} differs from the target Hilbert polynomial by a constant: to determine this constant we compare the 
value p{r) of the Hilbert polynomial p{t) in degree r with the cardinality of the order set .yV obtained by 
decreasing moves from .yV . 

Lemma 4.2. Let ,yY C V{n — l,r) be an order set and let jY C Vijijr) be the order set defined from ,jV by 
decreasing moves. Then, 

Proof. Each monomial G yV imposes the belonging to ^ (in addition to itself) of the monomials obtained 
from it applying the decreasing moves {e^ , 2e^, . . . , aie^}. □ 

There are three possibilities: 

• pi''') ~ \^\ < 0, ^ imposes too many monomials ouside the ideal, so the hyperplane section defined by 
^^^sat j^g^g discarded (there exist no Borel-fixed ideals corresponding to p{t) with such a hyperplane 
section); 

• p(r) — — 0, {SSy^"^ C K[xo, . . . , Xn] is one of the ideals sought; 

• ?'('') ~ \--^\ > 0, applying repeatedly Lemma 3.1 we determine the ideals we are looking for. 

Proposition 4.3 (Cf. [ Remark 5.3]). Let p{t) be an admissible Hilbert polynomial with Gotzmann number r 
and let ri be the Gotzmann number of Ap{t) . 
(i) Given a saturated Borel-fixed ideal J C K[a;i, . . . , a;„] such that K[a;i, . . . , J has Hilbert polynomial 
Ap(t), to pass from {(J • K\xq, . . . ,Xn])r} C 7^(n, r) to a Borel set corresponding to p{t), we need to 
remove at most r — ri monomials, 
(ii) We need to remove exactly r — ri monomials if we consider the lexicographic ideal L d K[a;i, . . . , a;„] 
corresponding to the polynomial Ap{t). 

Proof, (i) The minimal Hilbert polynomial having first difference equal to A.p{t) is T,{Ap){t) and the Gotzmann 
number of T,{Ap){t) coincides with the Gotzmann number of Ap{t). 

(ii) The Hilbert polynomial Y,{Ap){t) is admissible so there exists the saturated lexicographic ideal L with 
such Hilbert polynomial. By construction L is the 1-lifting of the ideal defining its hyperplane section, therefore 
it has no generators involving the variable xi and the ideal L n K[a:i, . . . , Xn], being generated by the same 
monomials, is still a lexicographic ideal. □ 

4.1. The pseudocode description. The recursive strategy naturally gives rise to a rooted tree where the 
nodes are all the Borel-fixed ideals met during the computation and the father of each node is the ideal of 
its hyperplane section (see Figure 1 for an example). This graph turns out to be a rooted tree because with 
a sequence of hyperplane section we obtain a unique ideal from any Borel-fixed ideal with the given Hilbert 
polynomial: 

• the ideal (1) C K[a:rf+i, . . . , x„] applying d + I sections, whenever the degree of the Hilbert polynomial 
d is smaller than n — 1; 

• the ideal {x'^) C K[a;„_i, a;„], where c = A'^p{r), ii d = n — 1. 



6 



P. LELLA 



The leaves are the Borel-fixed ideal we are looking for and by definition this tree will have height d + 2 if 
d < n — I and d + I if d = n — I. 

To compute the Borel-fixed ideals in K[a;o, . . . , x„] with Hilbert polynomial p{t), the algorithm described in [5] 
requires to compute the ideals in ]K[a;i, . . . , Xn] with Hilbert polynomial A.p{t) and so on, until the last recursive 
call that requires to compute the ideal representing the root of the tree. Now, starting from the root of the tree, 
the algorithm performs a Breadth-First-Search visit. This approach is not optimal, because it requires much 
space in the memory to store also all the ideal which are not leaves at maximal height and so which are necessary 
only temporarily. Moreover some of such ideals can be discarded immediately because too many conditions are 
imposed when embedded in a polynomial ring with one more variable. For this reason we prefer an algorithm that 
visits the tree with a Depth-First-Search approach. The function BorelFixedIdealsGenerator (Algorithm 
1) initializes the computation determining the root of the tree and starting the DFS visit. 



(1) 




{x^, X2X^) (x^, X2, X2X1, ^23;^) (x"^, x^X2, x^x-i, x^, x)^x'^){x'^, x^X2, x^x^ , x^x-^) 3:33:2, x^x^, X2){x'^, 3:33:2, 3:33:23:1, 3:33:2, 3:|) (3:3, 3:33:2, 3:2) 



Figure 1. The tree associated at the computation of the Borel-fixed ideals in lK[xo, 2:1,3:2, 2:3] 
with Hilbert polynomial p(t) — 5t ~ 2. 



1 


BorelFixedIdealsGenerator (IK[j:o, . . . , a;„ 






Input: IK[a;o, . . . , Xn], polynomial ring. 




Input: p{t), admissible Hilbert polynomial in P". 






Output: the set of all Borel-fixed ideals in IK[xo, . 


. , Xn] with Hilbert polynomial p{t). 




2 


d ^ degp{t); 






3 


r <— GOTZMANNNUMBER(p(i)); 






4 


if d — n — 1 then 






5 


c ^ A'^pit); 






6 


C7'(2,r); 




// K[Xn--i,Xn]. 


7 


return BORELlDEALs(K[a;o, . . . , x„],p(t), d, 






8 


else 






9 


^ ^ V(n - d,r)] 




II K[xd+i, . . . ,a;„]. 


10 


return BORELlDEALs(K[a;o, . . . ,x„],p(t),d 






11 


end if 







Algorithm 1: The procedure that initializes the computation detecting the root of the tree. 



The core of the algorithm is described by the function BorelIdeals (Algorithm 2). It takes a Borel set 
SS corresponding to a Borel-fixed ideal in K[a;fc, . . . , a;„] with Hilbert polynomial A''p{t). If fc = then 
corresponds to one of the searched ideal, so the function returns the saturation of the ideal generated by the 
monomials in Otherwise, the function embeds ^ in a poset with one more variable obtaining it computes 
how many monomials do not belong to ^ and how many monomials (q) have to be removed. Then it calls the 
function REMOVE to compute all the Borel set that can be obtained from ^ removing q monomials and finally 
for each of such new Borel set we have the recursive call with index k ~ 1. 

Also the function Remove (Algorithm 3) uses a recursive strategy, i.e. given a Borel set ^ and q monomials 
to remove, for each minimal element x" of it calls itself on the Borel set 3§ \ {x"} to which we have to 
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1 


BORELlDEALS(K[a;o, . . . , x„],p(t), fc, 




Input: ]K[a;o, . . . , a:„], polynomial ring. 




Input: p{t), admissible Hilbert polynomial in P". 




Input: k, integer s.t. ^ fc ^ degp{t). 




Input: W, Borel set in V{n - k, r) s.t. the ideal (l^)"""* in Kfxfc, . . . 


y ] Tl Q Q TTl 1 l^ll^T^" T\/~ll /^TTl 1 £1 1 /\ ^ Tif 

, lidiO ±liiUdlj JJUiy liUliiicLi Z_i jJ\L J . 


Output: the set of all saturated Borel-fixed ideals J in K[a;o, 


. . , Xfi] with Hilbert polynomial p{t) s.t. 




{Jr}{-^k) = ^■ 




2 


if fc = then 




3 


return {^f'^^; 




4 


else 




5 


^ { {{^Y^' ■ K[xfe_i, . . . , x„]) J ClV{n^k + l,r)- 




6 


q ^ A'^-VW - ; 




7 


if g ^ then 




8 


oneMoreVariable 4- Remove(^, q, 1); 




9 


ideals <— 0; 




10 


for all ^ G oneMoreVariable do 




11 


ideals ^ ideals U 






BorelIdeals(K[xo, . . . ,Xn],p{t),k 




12 


end for 




13 


return ideals; 




14 


else 




15 


return 0; 




16 


end if 




17 


end if 





Algorithm 2: The core of the recursive strategy to compute Borel-fixed ideals in K[x] with Hilbert polynomial 
Pit)- 



remove q — I monomials. Whenever we consider a Borel set ^ with (at least) two minimal elements x" and x'^ 
and we need to remove (at least) two monomials, the strategy just described would generate two times the same 
Borel set because (^ \ {a;"}) \ {a;^} = (^ \ {x"'}) \{x°'}. To avoid this repetition we add, as argument of the 
function REMOVE, a monomial x^ and we require that any monomial removed from ^ is greater than x^ w.r.t. 
a fixed term ordering (in Algorithm 3 we chose DegLex). So whenever the function is called by BorelIdeals 
we pass as argument the monomial 1 (any removal is admissible), whereas when Remove is called by itself, 
that is some removal has been already performed, we pass as argument the last monomial removed. In this way, 
assuming x" >DegLex x'' , the function Remove will only generate the Borel set \ {x''}) \ {a;"}. 

4.2. The implementation. The key of an efficient implementation is how we realize the Borel set. The 
guidelines are 

• slim structure, in order to take up less memory as possible (since the number of final objects can be 
huge); 

• quick implementation of the basic operations: 

(1) embedding of a Borel set in a poset with one more variable (lines 5-6 Algorithm 2); 

(2) computation of the minimal elements of a Borel set (line 6 Algorithm 3); 

(3) removal of a monomial from a Borel set (line 9 Algorithm 3); 

(4) computation of the saturation of the ideal generated by a Borel set (line 3 Algorithm 2). 

The most compact way to describe a Borel set is to store its minimal elements. This is efficient also to 
evaluate the belonging to the Borel set of any other monomial, because we can use the following proposition. 

Proposition 4.4 (['", Proposition 2.3]). Letx°'^x^ be two monomials m IK[a;o, . . . , x„]r. 

Tl 

x°'>BX^ ^ ^(aj -/3j) > 0, V i = 0,...,n. (7) 
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1 


Remove(^, 


Input: SS, a Borel set. 


Input: g, a non-negative integer. 


Input: x^, a monomial. 


Output: the set of all Borel sets obtained from SS removing in all the possible ways q monomials. The monomial 




is to avoid repetitions: it will be 1 when the function is called by BorelIdeals, whereas it will be the 




last monomial removed whenever the function is called by itself. 


2 


if g = then 


3 


return {.^}; 


4 


else 


5 


borelSets <— 0; 


6 


minimalMonomials ^ MinimalElements(^); 


7 


for all a;" G minimalMonomials do 


8 


if >DegLex then 


9 


borelSets borelSets U 




REMOVE(^\{a;"},q- 


10 


end if 


11 


end for 


12 


return borelSets; 


13 


end if 



Algorithm 3: The function to remove monomials from a Borel set. 



GOTZMANNNUMBER(p(t)) 

Input: a Hilbert polynomial. 

Output: the Gotzmann number of pit). 

MinimalElements(^) 
Input: a Borel set. 
Output: the set of minimal elements of 3§. 

Algorithm 4: Auxiliary functions 



Indeed if a monomial belongs to the Borel set, there exists a minimal element x" such that x^ >b x", 
whereas if x^ does not belong there exists a minimal element x^ such that x^ <b x^ . 

In this way, the computation of the set S of minimal elements is immediate. When a minimal monomial 
x" is removed, we computed all the monomials that can be obtained from x" with an increasing elementary 
move and for each of them we check if it belongs to the Borel set corresponding to S" \ {x"}: if not we add such 
monomial to the set of minimal elements. 

The description of the Borel set by means of its minimal elements turns out to be efficient also when we need to 
pass to a poset with one more variable. Indeed given the Borel set ^ C 7^(71 — fc, r) with minimal elements stored 
in the set M, by Lemma 4.1 we know that the Borel set SS = {{{My^'^ ■ K[xfe_i, . . . , x„])^} C V{n - fc + l,r) 
is described by the set of minimal elements 

M = {x^ . • . xlX\^xtL, I x«" . . . xlXX^xl'^ e M) . (8) 

A not trivial task is to compute the number of monomials that do not belong to the Borel set after the 
embedding in a poset with one more variable. In principle it consists in the computation of the Hilbert 
polynomial and in its evaluation at degree r. But the algebraic approach (used in [12]) slows down the algorithm 
because any software dedicated to the study of polynomial ideals computes the Hilbert polynomial with Grobner 
basis tools, which in this case are unnecessary. 



Given a Borel set ^ C Vin — k, r), let jV = . If we decompose jV according to the power of the smallest 
variable x^, i.e. 

r 

Zr = ]jZri, = {x"" eZ¥ s.t. ak=i] . (9) 

We can compute the cardinality of .jV — rewriting (6) of Lemma 4.2 as 

r r 

Hence we add to the structure describing a Borel set an array of r + 1 integers, such that the i-th index is equal 
to the number of monomials not belonging to the Borel set with the power of the smallest variable equal to i. 
It is easy also to deduce the array corresponding to a Borel set after the embedding in a poset with one more 
variable, indeed any € JV i implies the belonging to JV of « + 1 monomials 

_ an . . . ^"(0^,0 ri m: 

so that 



1^1 = Ei^. 



i=0 



Moreover when a minimal monomial x" is removed from a Borel set, we increase by one the index of the array 
corresponding to the power of the smallest variable in a:" . 

The last operation which has to be as quick as possible is the computation of the saturation of the ideal. We 
exploit the algebraic approach of [ ■ :] and we add to the structure describing the Borel set the list of monomials 
that generate the saturated ideal corresponding to the Borel set. If x" = x"" • • -x^*" is a minimal monomial of 
SS then x— = x"|a;^=i = x"" • ■ • x^j^-^ is a generator of the saturated ideal. As shown in [12], removing x" from 
^ implies that \ {x"})*''** is generated by the same generators of except for x— that is replaced by 

X— -Xm, ■•■ ,2^— -xi, where Xm = min X—. 



4.3. Experimental results. The following experimental results are obtained with an implementation of the 
algorithm realized with the guidelines showed above, coded in Java and ran on a MacBook Pro with a 2.4 GHz 
Intel Core 2 Duo processor. It can be tested by means of an applet available at 

www . personalweb . unito . it/paolo . lella/HSC/borelGenerator . html. 

In the following tables we consider Hilbert polynomials of degree 0,1,2 (i.e. points, curves and surfaces) and 
projective spaces of increasing dimension and we reported the elapsed time (in seconds) of the computation and 
the number of Borel-fixed ideals obtained. 
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Time (sec) 


n = 5 


n = 10 


n = 15 


n = 20 


Pit) = 5 


0.101 


0.025 


0.017 


0.021 


Pit) = 10 


0.062 


0.064 


0.119 


0.048 


Pit) = 15 


0.079 


0.225 


0.298 


0.401 


pit) = 20 


0.341 


1.595 


2.735 


3.870 


Pit) = 25 


2.094 


13.595 


24.497 


33.303 


Ideals 


n = 5 


n= 10 


n= 15 


n = 20 


Pit) = 5 


5 


5 


5 


5 


Pit) = 10 


42 


50 


50 


50 


Pit) = 15 


287 


417 


425 


425 


Pit) = 20 


1732 


3130 


3263 


3271 


Pit) = 25 


9501 


21616 


23158 


23291 



Table 1. Experiments with constant Hilbert polynomials. The Gotzmann number coincides 
with the number of points. 



Time (sec) 


n = 5 


n = 10 


n = 15 


n = 20 


5t + 1 (11) 


0.117 


0.159 


0.067 


0.0621 


5t + 7 (17) 


0.502 


1.480 


2.312 


3.290 


5i + 13 (23) 


10.513 


56.652 


91.456 


128.341 


8t - 6 (22) 


0.987 


2.623 


4.138 


5.852 


8< - 3 (25) 


3.008 


14.128 


22.960 


32.300 


8t (28) 
Ideals 


12.960 
n = 5 


72.053 
n = 10 


273.719 
n = 15 


238.856 
n = 20 


5t+l (11) 


89 


98 


98 


98 


5t + 7 (17) 


3028 


4560 


4587 


4587 


5t + 13 (23) 


58124 


123689 


126962 


127030 


8t-6 (22) 


4171 


6741 


6837 


6838 


8t-3 (25) 


17334 


32073 


32848 


32868 


8t (28) 


68291 


144660 


149777 


149976 



Table 2. Experiments with Hilbert polynomials of curves. The Gotzmann number is reported 
in brackets. 
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Time 1 


sec) 


n = 5 


n = 10 


n = 15 


71 = 20 




+ 8*- 


46 (16) 


0.312 


0.189 


0.304 


0.516 


2*2 


+ 8*- 


42 (20) 


0.103 


0.338 


0.558 


0.883 


2*2 


+ 8*- 


38 (24) 


0.741 


3.167 


5.237 


7.055 


4*2 


- 12*4 


- 10 (20) 


0.147 


0.280 


0.377 


0.561 


4*2 


- 12*4 


- 14 (24) 


0.953 


3.909 


6.007 


8.588 


4*2 


- 12*4 


- 18 (28) 


9.066 


50.071 


82.592 


112.237 




Ideals 


Tl 5 


77,-10 


Tl 15 


71 = 20 


2*2 


+ 8*- 


46 (16) 


834 


38 


38 


38 


2*2 


+ 8*- 


42 (20) 


481 


670 


671 


671 


2*2 


+ 8*- 


38 (24) 


4774 


8393 


8476 


8476 


4*2 


- 12*4 


' 10 (20) 


631 


856 


857 


857 


4*2 


- 12*4 


- 14 (24) 


6394 


10986 


11082 


11082 


4*2 


- 12*4 


- 18 (28) 


51527 


112852 


115295 


115332 



Table 3. Experiments with Hilbert polynomials of degree 2. The Gotzmann number is 
reported in brackets. 
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